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Abstract 

Neutrinoless double beta decay data together with information on the absolute neutrino masses 
obtained from the future KATRIN experiment and/or astrophysical measurements give a chance 
to find CP violation in the lepton sector with Major ana neutrinos. We derive and discuss necessary 
conditions which make discovery of such CP violation possible for the future neutrino oscillation 
and mass measurements data. 

PACS numbers: 14.60.Pq, 26.65. +t, 95.85.Ry 
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I. INTRODUCTION 



Information on CP violation in the lepton sector is very important for building the future 
theories which go beyond the Standard Model As CP violation is probably predomi- 
nantly connected with lepton masses and observed neutrinos are very light, an experimental 
measurement of the effect is a serious challenge. For three Dirac neutrinos there is one CP 
violating phase (6) and additional two phases (ai,tt2) exist for Majorana neutrinos. The 
charged current state (ua) is related to mass states (z/j) by an unitary transformation 
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C12C13 
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V 

Cij and Sij are cosines and sines of the 9ij {ij = 12, 13, 23) angles. The second matrix in (j21) 
appears only for Majorana neutrinos. 

It is commonly believed that CP violation owing to Dirac phase S will be discovered in the 
future superbeam or neutrino factory experiments where oscillations of neutrinos and 

antineutrinos will be observed. From the parametrization of the mixing matrix (Eq|2) we 
can see that sin ^13 and e^*'' always appear in a combination. So, any CP breaking effect for 
Dirac neutrinos will be proportional to sin ^13 sin 5 and disappear for sin ^13 —>■ 0. From the 



present fits it follows that this mixing angle is small (sin^ ^13 < 0.05 for 99.7% C.L. 
and the assumption that ^13 = agrees with the data equally well. Such a tendency, if 
outlives in the future when more precise data will be available, effects in a very small or 
vanishing CP breaking signal. It was shown that for 5 = ±| effects of CP violation will be 
seen in the future experiments if sin^ 6'i3 is not smaller than 10~^ ||7|. 

If neutrinos are Majorana particles, in addition to the phase S, two other phases can 
also be responsible for CP symmetry breaking. Many different processes are, in principle, 
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sensitive to these Majorana phases and can generate both CP-even and CP-odd effects j8[|9|. 
Admittedly, most of them are much beyond an observable limit. The only experiment which 
could provide evidence for Majorana phases is the search for neutrinoless double beta decay 
Such a possibility has been discussed many times jiol . 111 . 



18 



19 



20 



21 



22 



15 



17 



23( 1 but, to our knownledge, detailed conditions concerning the future 



experimental results and their necessa ry p recision to discover CP violation has not been 
discussed. The exception is where Ref. j2J| authors consider the future anticipated precision 
of all relevant neutrino experiments, and formulate a very pessimistic 'no-go' conclusion. 
They state that even under a very optimistic assumption about sensitivity of the future 
experiments it will be impossible to detect neutrino CP violation in the {Pf3)ou decay. We 
agree with such a statement, but we would like to go a step further. We propose a set 
of conditions for neutrino masses and mixing angles (best fit values = b.f.v) altogether 
with conditions on experimental and theoretical precision for their determination, such that 
discovery of CP violation arising from Majorana phases in the {PP)ou decay will be possible. 
We formulate sufficient conditions when CP violation could be observed. We should mention 
that our conditions are completely general. Contrary to Ref. [24] we do not assume from the 
beginning that the 6*13 angle vanishes. Similar consideration has been done in 25, 26, [27 1. 
Here we concentrate on the degenerate neutrino mass spectrum where CP violation has a 
clear meaning. We investigate in more details the problem of theoretical determination of 
the nuclear matrix elements, mechanism responsible for [(3(3)o„ and the future experimental 
error of nuclei decay lifetime. 

We found that under a very optimistic assumption on the sensitivity of the future exper- 



iments considered in Ref. 



2J|, independently of measured b.f.v., it is really impossible to 



detect CP violation. However, such a possibility is 'just around the corner'. A little better 
precision will give a chance to make a decisive statement about CP Majorana breaking. 
Even if required precision for today is estimated to be very optimistic value, we hope that 
the problem of lepton CP violation is so important that it is worth to have it in mind. 

Other important result of our investigation concerns the ^13 and du mixing angles. In 
contrary to neutrino oscillation experiments, smaller ^13 angle gives better prospect of CP 
symmetry breaking measurement. Similar situation take place for solar mixing angle 612- 
The maximal 612, sin 2^12 —>■ 1, is the best scenario for CP symmetry breaking discovery, 
contrary to the prospects of finding the neutrino mass bound from {[3(3)qu decay j28!|, where 
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012 — ^ f brings such a possibility to ruin. 

In the next chapter we discuss how CP symmetry braking could be determined from 
neutrinoless double beta decay. Then, in Chapter III, we describe the present situation 
and we predict how precisely all parameters (oscillation mixing angles, effective mass {rriy) 
measured in (/3/5)oj/ and mp measured in e.g. tritium beta decay) should be determined in 
order to discover CP symmetry breaking. Two kinds of presentations are given. The first 
one, very visual, where correlations between errors are not included. And the second, where 
more sophisticated analysis which answer at which value of confidence level the probes of 
CP violation could be carried on. Finally, Chapter IV contains our conclusions. 



II. CP SYMMETRY BREAKING AND THE {(3(3)o„ DECAY. 

The neutrinoless double beta decay {f3P)ou of nuclei measures the effective neutrino mass 

3 
i=l 

where 02 = «2 — cti and 03 = — 5 — «i. 

As we will see, the possible precision of the future experiments will give a chance to look 
for CP violation only for higher neutrino masses (mi > 0.1 eV). For this case the mass 
spectrum starts to be degenerated and we will consider only such a spectrum. Then the 
effective neutrino mass rrif^ measured in tritium beta decay, independently of its definition 
is just equal to neutrino masses 

= 

For Majorana neutrinos CP symmetry holds if and S take one of the values 0, ±|, ±7r. 
Then from Eq. Q four conserving CP values of (m^) are obtained: 

{'mu)(i) = rrip, 

{mu){2) = m^cos26'i3, 

{'mu){3) = (cos^0i3| cos 26*12 1 + sin^^i3) , (5) 

("T'i/)(4) = (cos^ 0i3| cos2^i2| -sin^^i3) . 



i=l 



Elf. 



rrii 



mi ~ 1712 ~ m^. 



(4) 



4 



In all cases, the relation between {niv) and mp is linear with different slopes {i = 1,2,3,4) 

{mu)(i) = Ciirtp. (6) 

First we would like to present very visual method of finding a region of parameters where CP 
violation can be probe. Lately we will present a result with the correct statistical analysis. 
Let us assume that 9ij mixing angles are known with definite precision, 

sin^ % G ((sin^ Oij)min, (sin^ Oij)max) (7) 

with central value 

(sin 6ij) 

best fit • (8) 

For each q = 2, 3, 4) we can calculate the maximal and minimal values 

cr" = (cos 2^13)^ax, 
cr" = (cos 2^13)™.., 

C^«^ = (cOS^^i3)^ax(cOs2^i2)max + {sin^ 9 13) max, 

CT"" = (COS^ ^13)mm(cOs2^i2)™„ + (slu^ ^13)mm, (9) 
^rnax ^ [cOS^ di'i)max{cOS 26 12) max - (siu^ ^13)mm, 

C™" = (COS^ ^i3)mm(cOs2^12)mm ~ (siu^ ^13)max • 

We can see that localization of the (mj,)(j) lines is fully determined by the oscillation pa- 
rameters, namely ^13 and 612 angles. 

Let us now assume that in the future experiments m^j and {rriy) masses are determined 
with precision Am^ and A{m,^): 

{m^)exp ± A(mi,), (10) 
imp)exp ± Am/3. (11) 

Then localization of the rectangle R = (Am^j, A(mi^)) between the lines ci = 1 and c™" 
(see Fig^ decides about CP symmetry breaking. If R crosses the error region between the 
(c™", c^"^) lines z = 2, 3, 4, we do not know anything about CP symmetry. But, in opposite, 
if R is located outside the Cj error region then there is indication that CP symmetry is broken 
as at least one of the angles 6, ai, a2 is not equal to its CP conserving value. 
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FIG. 1: A localization of the R = {Ami3, A{m^)) rectangle between c™" and c™*^^ lines which 
indicates that CP symmetry is broken. 

Possible localization of the present and prospective (m,^)j = Cimp lines is presented in 
Fig. 121 and El respectively. We can see that localization of R between c^"^' and c^*" lines is 
only interesting for CP violation search. If the rectangle R with Am/3 and A{m^) sides is 
fully located between two lines with the c™""^ and c™*" slopes then CP symmetry is broken 
(see Fig^). So first conditions for detecting CP violation are: 

Am/3 < L, A(m^) < K. (12) 

L and K can be find in a easy way 

K = {m/3)A - iAm/3)B, (13) 

and 

L = {m,)C - A{m,)D, (14) 
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FIG. 2: A localization of the (c™" c^"^) regions for the present ^13 and 612 angles precision. To see 
CP violation a precision of mp and {m^) measurements should be very good. For smaller (and 
{my)) a region where CP violation can be search for is smaller, so a precision of their measurements 
should be even better. 

where 



A mm max 

/I — 62 — C3 , 



B 



„mm I „max 

-2 "T "-a 



C- 
D 



A 



„min„max ' 
^2 



(15) 



B 



„max 
-3 



for any and {rriy) values inside the two lines c^™ and Cg 

If conditions (Eq. [T21) are satisfied for some central values {mp)exp and {my)exp determined 
from experiments (and theory) then there are further two possibilities. The rectangle R 
located in the point ((m^)ea;p, {innu)exp) can 
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FIG. 3: The CP conserving regions (hatched areas) which follow from the future neutrino oscillation 
experiments. We assume that central values of ^13 and 612 are in agreement with present data 
but their error estimation is suppose to be much better (sin^ 612 = 0.28 it 0.01 and sin^ ^13 = 
0.005 lb 0.0001). The region between c™" and c^°-^ lines is larger, giving more space for the 
rectangle (Am/3, ^{m^)) (see text for more details). 

1. be fully inside two bounding lines c^**^ and c^"'^, or 

2. be located partly on the first or the second line. 

In the first case we can conclude that CP symmetry is broken, in the second the problem is 
unresolved. The first condition is satisfied if: 

cT' {{mpU + ^) < [{m.)e., - ^) , (16) 

and 

((-.)e.p + ^) < {{m,U - ^) • (17) 
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The inequalities given by Eqs. ffT^ . (fT^ and (fT7|) form the set of necessary conditions for 
CP symmetry breaking. Of course we are not able to prove in this way that CP symmetry 
holds. 

Let us parameterize 

A{my) =2x{mu), Amf3 = 2ymf3, (18) 

where 2x is the relative error which measures the theoretical nuclear matrix elements uncer- 
tainty and experimental decay lifetime of the (m^) matrix element. Similarly, 2y measures 
the relative error of the effective mass e.g. from tritium beta decay. As both K and L (in 
Eq. (fT^ and (HH)) must be larger than zero, we have two consistency conditions. Both x 
and y must satisfy the same inequality 

1 — cos 2^12 min — 3 siu^ ^13 max + siu^ ^13 min COS 2^i2 min 

•^^ y — ~, I 7^ • 2 n • 2 n ^W) ■ V / 

i + COS ZU12 min ~ t/l3 max "~ Sm f^lS min COS ZU12 min 

These inequalities impose sharp conditions concerning a precision of the and {rrii,) de- 
termination. As the r.h.s. of Eq. is a decreasing function of sin^ ^13 and cos2^^i2, the 
best circumstances arise for sin^ ^13 and sin^ ^12 — ^ |- In this case lines c^*" 1 and 
^max _^ g gj^g ^^^q largest region for localization of {mu)f.xp and {mfj)f.xp where symmetry 
is broken. As we know, the condition ^13 ruins the Dirac 5 phase determination in 
oscillation experiments. We can see that both methods, {f^P)ou decay and long baseline 
experiment which could detect 5 are complementary for detecting CP violation j^. Also 
the other condition, the large solar mixing angle [612 —>■ j) is not favourable for Majorana 
mass determination from the (/3/3)oi/ decay [7-18]. 



The case sin^ 6*13 = has been considered in Ref. 



24. Then Eq.^ gives 



X < tan^ 9i2 (20) 



which is exactly the condition given by Eq. (14) in Ref. j 

From Eqs. and (fT^ for given relative errors x and Am^ we can also find the lower 
limit for the m/j and {rriiy) effective masses for which measurements are still possible 

> (21) 



and 

Z 

'^^^ ~ A ' C-2xD 



m,>^(B + -^^]. (22) 



Now, using present and expected in the future precision of the neutrino oscillation data 
we can estimate how well and (m^) should be determined to discover CP symmetry 
breaking. 



III. NUMERICAL RESULTS 



Using presently determined 9i2 and ^13 mixing angles 0, H Q, S Q (with 3cr preci- 
sion) 



from Eq. ()19|) we obtain: 



0.22^ sm^ei2 ^0.37, 
^ sin^ 2^13 ^ 0.048, 



X < 0.2. 



(23) 



(24) 



It will be a serious challenge to get such a precision. Let us check it for the isotope of 
Germanium ^^Ge where evidence for the (/3/3)ov decay is claimed to have been obtained |37|. 
If we assume that only one standard mechanism, the exchange of Majorana neutrinos with 
masses mj, is responsible for the (/3/3)oiy decay, the effective mass {rriy) is calculated from 
decay rate T{^^Ge) IsS]: 

T-\'^Ge) = G|M|2(m,)2, (25) 

where G is accurately calculable phase space integral and M is calculated Nuclear Matrix 
Element (NME). Unfortunately, this calculation is a complicated job, and different methods 
of calculation give different results. For the isotope ^^Ge the results differ by one order of 
magnitude. If we parameterize 



T(^6G'e) = h X 10^4 y, 
then for {mj) = 1 eV eleven different results have been obtained 

3.15 
17.7 



(26) 



b=1.7{=bmin) 2.16 41] I, 2.3 |41|L, 2.33 

3.6 H, 4.06 45|, 8.9513. 14.0 

However, we would like to stress, that method used in Ref. 



3.2 



421 



3,0, 



41 



44 



45 



46 



are completely independent, different nuclear models are used, and generally models 
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FIG. 4: The lines for the full uncertainty of {rriy), x = ^(!^"\^ as a function of both, theoretical 
uncertainty in NME calculations a = (see text) and experimental relative error (xt) for the 
decay lifetime of ^^Ge. 

are not calibrated against nuclear properties. If we assume that relative experimental error 
for TiJ^Ge) measurements is defined by 2xt 



2x1 



then the full relative uncertainty of (m^) {2x = ^^^) is given by {xt < I] 

+ XT)a - y/1 - Xt 



(27) 



X = x{a, Xt) 



(28) 



where a 



a/(1 + XT)a + y/1 - Xt 
^oiasL^ In Fig. 4 we can see the relation between the NME precision (a) and 
the expected uncertainties for effective neutrino mass {rriy) = ^^^^^^ various future 
experimental errors of the decay lifetime of ^^Ge. 



We can see, that taking seriously present discrepancy in the NME determination (a ^ 10) 
we obtain x ~ 0.52, much larger than necessary (see Eq. ()24|) ). The new calculation of NME 
where the observed {(3(3)2y decay has been used to fix relevant parameters, has shown 



11 



the great stability of the final results. For the ^^Ge two methods of calculation, RQRPA 
and QRPA (see j42] for more details) have given almost the same results, and then: 

a « 1.4. (29) 

With such a precision of the NME determination we obtain {xt ~ 0.3) 

X « 0.24, (30) 

still above the present necessary precision (see Eq. (j24|) ). but within reach of the future 
oscillation experiments. 

We should also mention the other uncertainty in the (m^) determination - the possible 
different physical mechanism for the {Pf3)ou decay. If the lepton number is violated at 
TeV scale we can expect the other processes which give equally strong, as light Majorana 
neutrinos exchange, contributions to {(3l3)ou- Then the relation between decay lifetime and 
{rriv) is not given by Eq. ()25j) . To answer a question at which scale lepton number is violated, 
information from higher energy colliders (e.g. LHC) and other lepton processes is necessary. 
In Ref. js^ it was shown that a study of two lepton flavor violating processes /i — e 
conversion and ^ e + 7 decay will give important insight to the mechanism of the {fip)^^ 
decay. 

From Eqs. ()2Hl and ()22|1 we can flnd conditions for m^g and {niv) effective masses for which 
CP symmetry breaking could be seen (see Fig(j2I)). For example, ii x 0.15 with the present 
3cr precision of mixing angles Eq. lj^Hj) and for Am^ = 0.03, 0.02, 0.015 eV, the CP symmetry 
breaking is testable for {rriy) > 0.24, 0.16, 0.12 eV and m/j > 0.32, 0.21, 0.16 eV, respectively. 
There is some chance that in the future experiments such Am^ precision can be reached, 
but relative error for {m,y), a; ~ 0.15 is much beyond the present possibilities. 

From Eq. (PT|) for a given central value of {rriiy), we can flnd relation between the x and 
Am/3 required for probing the CP symmetry breaking. Let us assume that a value of {rriy) 
is really in an interval given by the Heidelberg group 37 1 

{m,),,p ^ (0.1 - 0.9) eV. (31) 

If {m^)exp ~ 0.1 (0.9) eV, Anifs should be smaller than 0.002,0.013,0.026 (0.014,0.11,0.24) 
eV for X ~ 0.19, 0.15 and 0.1, respectively, with the central value m/3 ~ 0.13 (1.2) eV. 

More careful analysis, taking into account the present precision of the mixing angle de 
termination 



i33 | can give a region in the {{m^),mp) plane where CP violation can be probe 
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FIG. 5: The regions of relative error of rrii, versus m^, where CP is violated with confidence level 
(CL) equal to: 90%, 95%, 99% and 99.97% {3a). To find CP symmetry breaking, {rriu) and 
should be determined with extremly difficult to reach precision. 

with various CL. The regions of relative errors and for which CP violation could 
be seen are presented in Fig. 5. We see that even for 90% CL the x parameter should 
be smaller than x < 0.15, so it is completely out of reach with present experimental and 
theoretical possibilities. 

How a better determination of the Ou and ^13 mixing angles affects the x and Am^ 
uncertainties? Let us assume that during next years the precision of experiments will be 
strongly improved. Let us also assume that the best values of mixing angles will not change 
but only precision will be much better: 

1. The KamLAND and Borexino experiments determine the solar mixing angle with 
precision sin^ 612 ~ 0.28 ± 0.01 



2. The IHF-Kamioka neutrino experiment or the future neutrino factories [53 will mea- 
sure the with the precision A^i3 = 0.01 (so sin^ = 0.005 ± 0.0001). 
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End assume finally that: 

3. Weak lensing of galaxies by large scale structure together with CMB data measure the 
sum of neutrino masses Yl = ""^i + + 1^3 to an uncertainty of 0.04 eV. So we can 
expect that each individual mass is known with the precision Am/3 = 0.015 eV js^. 

Now from Eq. (jl9p we get the required precision of Arrip and A{rnJ) 

X, y < 0.36. (32) 

In Fig. 4 we present for this value of x a necessary precision of NME for different relative 
errors of the T(^^Ge) measurements. If the last estimation of NME is confirmed (a » 1.4) 
and the decay lifetime of ''^Ge is found with xt < 0.5 then necessary precision of (m,^) will 
be obtained. Such a scenario is not only a pure fantasy. More precise estimation will give a 
region of {m^) and where a probe of CP violation could be possible (Fig. 6). We have 
assumed the same relative uncertainties for (m,^) and m/? {x = y). For x = y = 0.07 there is 
no region where CP could be find with CL > 99%. This region appears if x = y are smaller. 
In Fig. 7 a region of x and y relative errors is presented for a given level of CL. We can see 
that if we want to probe CP violation with CL ^ 90% x must be smaller than x < 0.22 for 
very well determined nip {y 0) and vice versa, y < 2 for x 0. Correlations between 
quantities give more stringent requirements for relative errors (see Eq. ()32|) ). We can see 
from Fig. 4 that to get x ~ 0.1, parameter a must be smaller than 1.3 and xt better than 
10%. Knownledge of NME on a 30% level has been postulated recently 



IV. CONCLUSIONS 



From presented estimations it follows that measurement for CP violation for Majorana 
neutrinos in neutrinoless double beta decay could be possible for almost degenerate spectrum 
of their masses (m/j > 0.1 eV). However, several conditions should be satisfied: 

1. oscillation mixing angles should be measured with better precision e.g. A(sin^i3_« 
0.01 and A(sin6'i2 ~ 0.1) which are within the future experimental range (see e.g.jSl], 

2. absolute neutrino masses ruf^ should be measured with precision Amp ^ 0.02 eV with 
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FIG. 6: Regions in the (mj^) <^ rap plane where CP symmetry is broken with various CL for equal 
relative errors of {my) and (x=y). 

.he central value in the range mp > 0.15 eV, which is also not a fully fantastic dream 

3. neutrinoless double beta decay is discovered and the decay lifetime T is measured 
with precision better than 10%. It is difficult to say at the moment anything about 
the future precision of T. If we give a credit to the last Heidelberg group news about 
{PP)ou decay of ^^Ge, then the error of T is much higher. They derived from the full 
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FIG. 7: The regions of relative errors of {m^) versus where CP is violated with CL equal to: 
90%, 95%, 99% and 99.97% {3a). Uy = 0.05 then to determine CP violation, at 99% CL we have 
to know the effective Majorana mass with precision x = 0.09. 



data taken until May 2003 that [3 

TC^Ge) = (0.69 - 4.18) x 10^^ y, 

with the best value TiJ^Ge) = 1.19 x 10^^ y, so the relative error xt 
get Xt < 0.1 will be probably a very difficult task. 



(33) 



^ ~ 2.9. To 



4. nuclear matrix elements of decaying isotopes are calculated with much better precision. 
Future uncertainties for a = ¥^ should be smaller than a < 1.3. During the last 
years some improvement in NME calculation has been obtained. The last result where 
a ~ 1.4 has been presented is a very good step forward j49^. The model of NME 



calculation can also be tested via comparison o: 
(or more) nuclei with experimental data 0,0, 
(/3/?)oiy decay of several nuclei is observed. 



t 



le results of calculation for three 



57| . This test can be acomplished if 



5. there should be independent information about a full mechanism of the {Pf3)oi, de- 



Id 



cay. We should know that two electrons are produced by two W-bosons and Majorana 
neutrino exchange virtual process. Any other mechanism should give negligible contri- 
bution to the neutrinoless electrons production. The future LHC data and observation 
of other lepton violating processes give some chance to clarify this issue. 5ol |. 
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